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NATIONAL ADVISORY COMMITTEE FOE AERONAUTICS 

TECHNICAL NOTE NO. 1269 

^^ETHOD FOR CALCULATING '.VENG CHARACTERISTICS 
BY LH’TING-LINE THEORY USING NONLINEAR 
SECTION LIFT DATA 

By James C . Si veils and Robert H • Neely 
SUNflARY 


A method is presented for calculating wing characteristics by 
lifting— line theory using nonlinear section lift data. Material 
from various sources is combined with some original work into the 
single complete method described, Multhopp'e systems of multipliers 
are employed to obtain t)ie induced angle of attack directly from 
the spanwlse lift distribution. Equations are developed for 
obtaining these multipliers for any even number of spenwise stations, 
and values are tabulated for ten stations along the semlspan for 
asymmetrical, sjTumetrical, and antisymmetrical lift distributions. 

In order to minim.ize the computing time and to illustrate the 
procedures involved, simplified computing forms containing detailed 
examples ere given for symmetrical lift distributions. Similar 
forms for as;anmetrical and antisymmetrical lift distributions, 
although not shown, can be readily constructed in the same manner 
as those given. The adaptation of the method for use v;ith linear 
section lift data is also Illustrated. This adaptation has been 
found to require less computing time than most existing methods. 

The wing characteristics calculated from general nonlinear 
section lift data have been found to agree much closer with 
experimental data in the region of maximum lift coefficient than 
those calculated on the assumption of linear section lift curves. 

The calculations are subject to the limitations of lifting-line theory 
and should not be expected to give accurate results for wings of 
low aspect ratio and large amoimts of sweep. 


INTRODUCTION 


The lifting— line theory is the best known and most readily 
applied theory for obtaining the spanwlse lift distribution of 
a wing and the subsequent determination of the aerodynaiilc character- 
istics of the wing from two-dimensional airfoil data. The character- 
istics BO determined are in fairly close ag 2 -eement with experimental 
results for wings with small amounts of sweep and with moderate 
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to high values of aspect ratio; for this reason, this theory has 
served as the hacis for a lai’ge part of present aeronautical 
hnowledgo . 


Tie hypothesis upon which tho theory is based is that a 
lifting wing can be replaced by a lifting line and t'lat the 
incremental vortices shed along the span trail behind the wing in 
straight lines in the direction of the free-stream velocity. The 
strength of these trailing vortices is oroportlonal to the rate 
of change of the lift along the span. The trailing vertices induce 
a velocity normal to the direction of the free— stream velocity and 
to the lifting line. The effective angle of attack of each section 
of the wing is therefore different from the geometric angle of 
attack by the amount of the engle (called the induced angle of 
attack) whose tangent is the ratio of the value of the Induced 
velocity at the lifting line to the value of the free— stream 
velocity. The effective angle of attack is thus related to the 
lift distribution through the induced angle of attaclc. In addition, 
the effective a.ngle of attack is related to the section lift 
coefficient according to two-dimensional data for the airfoil 
sections incorporated in the wing. Both relationships must be 
simultaneously satisf 1 ed in the calculation of the lift distribution 
of the wing. 

If the section lift curves are linear*, these relationships 
may be expressed by a single equation which can be solved analytically. 
In general, however, the section lift curves are not llneai*, 
particularly at high angles of attack, and analytical solutions 
are not feasible. The method of calculating the spanwise lift 
distribution using nonlinear section lift data thus becomes one of 
making successive approximations of the lift distribution until 
one is found that simultaneously satisfies the aforementioned 
relationships. 

Such a method has heen used by Wieaelsberger (reference 1) 
for the region of maximum lift coefficient and by Boshar (reference 2) 
for high— Bubsonic speeds. Beth of these writers used Tani’s 
system of miltipliers for obtaining the induced angle of attack 
at five stations along the semlsprn of the wing (reference 3 ) • 

Tanl, however, considered only the case of wings with symmetrical 
lift distributions. Multhopp (reference 4), using a somewbjit 
different aatheimitico.l treatment froBi that which Tanl used, derived 
systoms of multipliers for syrametrical, anti syraaetri cal, and 
asymmetrical lift distributions for four, eight, and sixteen 
stations along the semispan. Multhopp‘s derivation, in slightly 
different form and ncmenclature, is presented herein and tables 
are given for the multipliers for ten stations along the semispan 
(the usual number of stations considered in i£ifin.y reports in the 
United States). 
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For eymmetrical distributions of wing chord and angle of 
attack^ the multipliers for symmetrical lift distributions may be 
used with nonlinear or linear section lift curves. For asymmetrical 
distributions of angle of attack^ the multipliers for asymmetrical 
lift distributions must be used if nonlinear section lift curves 
are used. If an asymmetrical distribution of angle of attack can be 
broken up into a symmetrical and an anti sr'mmetri cal distribution, 
the anti symmetrical part may bo treated separately if the section 
lift curves can be assumed to be linear. 

The purpose of the present paper is to combine the ccntrlbutloiB 
of Multhopp and several other writers, together with some original 
work, into a single complete method of calculating the lift 
distributions and force and moment characteristics of wings, using 
nonlinear section lift data. Simplified computing forms are given 
for the calculation of symmetrical lift distributions and their use 
is illustrated by a detailed example. The adaptation of the method 
for use wi.th linear section lift data is also illustrated. No 
forms are given for asymmetrical or anti symmetrical lift distributions 
inasmuch as such forms would be very similar to those given. 


SYMBOLS 


S 

b 

c 

^s 

^t 

c* 

c» 

A 

X 

y 

z 


wing ai"ea 
wing span 

chord at any section 
root chord 
tip chord 

mean geometric chord (s/b) 

po/2 2 

cLy] 

\ s I 

aspect ratio (b^/s) 

coordinate parallel to root chord 

coordinate perpendicular to plane of symmetry 

coordinate perpendicular to root chord and parallel to 
plane of s^OTnetry 


mean aerodynamic chord ( — 

I S 
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free— etreem dynamic presBUr© 

■Reynolds nuaibor 

mass density 

free— stream velocity 

coefficient of viscoaity 

wing lift coefficient (l/qS) 

section lift coefficient (l/qc) 

wing lift 

section lift 

wing drag coefficient (C/qS) 
wing profile-drag coefficient 
wing induced-drag coefficient 
section prof ile— <drag coefficient 
section induced-drag coefficient 
wing drag 

wing pitching-HHoment coefficient (M/qSc’) 

section pitching— dOTnent coefficient about section quarter- 
chord point 

wing pitching moment 

wing rolling-moment coefficient (L*/qS) 
wing rolling moinent 

wing induced— yawing-moment coefficient 

wing prcfile-yawlng^Roment coefficient 

anglo O''"' attack of any section along the span referred 
to its chord l;'no 
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angle of nttaolc of root section referred to Its chord 
lino 

angle of attack of root section referred to its zero 
lift line 

section induced angle of attack 

effective angle of attack of anj’' section 

section angle of attack for tvo-dinensional airfoils 


ngle of zero lift of any section 




angle of zero lift of root eoct:‘ on 


wing angle of atcack for .eero lift 

geometric angle of tvist of any eectlon along the epan 
(negative if washout) 

aerc'iyra'^mic arglo of tvrist of any section along the span 
(negative if washout) 

gocmet'vic angle of tvist of tip section 

aercfynarxilc angle of twist of tip section 

wing lift— curve slope, per degree 

section lift— c’jrve slope, per degree 
I T''VO'-d.J.mensional lift— curve slopo\ 

\ Edge-velocity factor j 


COB Q coordinate (Sj/t) 


coefficients in trigonoiDetric series 

TiiuI.tlplior for Induced angle of attack (asymmetrical 
distributions) 

multiplier for induced angle of attack (symmetrical 
distributions) 


multi rlier for induced angle of attack (antisymmetrlcal 
distributions) 
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Tj multiplier for lift, drag, and pitching-moment coefficients 

(aa^mmetrical distritutions) 

T) mulbiplier for lift, drag, and pitching— moment coefficients 

(sj^rametrical dlstrihu-ciona) 

0 ^ multiplier for rolling- and yawing-moment coefficients 

( as’mmctr i cal di ctritutl ons ) 

0 mult?pl-^er for rolling-moment coefficient (anti sj^rnmetri cal 

di strlbutiona ) 

E edg.o-volocity factor 

\ span / 

Subscripts 

max iLaxiErm va] ue 

al Talue for additional lift (Cj^ = l) 

b value for basic lift (Cg = O) 

/a„ ^ value for conctant value of a„ 
val'ue for given value of 


TERIOIRETIC/iL DEVELOPMENT OE METHOD 


Lift Distribution 


The methods of Tanl (reference 3 ) and Multhopp (reference 4) 
for determining tlio induced angle of attack are fundamental!?;- 
the same, differing only in the mathematical treatment. The 
method presented herein is essentially the same as that given 
by Multhopp, In tlie following deidvation t]->.e spanwise lift 
distribution is e:q>ressed as the trigonometric series 

-J- - Bin no (1) 

as in reference 9, vrhere 9 is defined by the relation cos 0 = , 

It ma,v be noted that each coefficient ,1^^, as used lierein, is equal 
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to fo'ir tJmgs the correepcndlng coefficient in reference 5, The 
inluced angle of attack (in defprees) at a point y, on the l:'fting 
line io 





pW'2 /111) 

13 Q ^ 

h J 

7T 8jT 

' 

u 

-b/2 


yi ■ ■ y 


dy 


(2) 


This integral (in different nomenc J atm-e ) was given hy Prandtl in 
reference 6, If equation (l) is suostituted into equation (2) and 
the variahle is chonged frori y to Q, the induced angle of attack 
at the general point 9 teooriyjs, according to reference 


a = - 

4it sin 0 --- 


nAj^ sin nS 


(3) 


The prchlom of obtaining the induced angle of attack is thus reduced 
to cne cf d otcnulning the coefficients of the trigonometric series. 


Tile lift di str:' ’mti on (equation (l)) may Te approximated ly 
a fin: ue tr ■gciiometric aeries of r — 1 terms whores for suhoequent 

c c 

usa^e, r \s asc'jiricd to even. The values of — at the ccinally 

b 


epaceJ. po:?ntB 0 


mit 

r 


in the ra^^yo 0 < 0 < jt aue expi^essed as 



(^) 


c « c 

vhere n = 1, 2^ r — 1. Conversely^ if the valuos of — ^ 

b 

are kno\ni at each point the coefficients of the finite sei*ies 

may bo founl l:y hamonic analysis as 



(?) 



8 


KACA TN No . 1269 


Tf equetion (5) is Bu’bstitu'ted In equation (3)j a doutle 
summation is obtained for the induced angle of attack as 


0^(9) = 


l8o_ 

Bin d 



n sin ne 


r— 1 

2 V 

^ A- 

El=l 


C-i^) e,n 

\ t /m 


5? 

r 



If fhe ii? 4 uc©d ©iigle of attack ic to bo detormlned at tho sam© 
points 6 rt ‘Kbich the load distribution is known^ that is^ at the 
points 0 = then 


“Iv = 


180 


r-1 


(k — m)it (k + in)rt 

cos n cos n 1 


r-1 



mp:l 


( 6 ) 


vhero 


P mk ~ 


r-1 

, - ^ n 

4Jtr Bin ^ 

'' n=l 


180 


cos 


n _ 


COB n 


(k + m)jt 
r 


( 7 ) 


It can be shovm that, if cob ^ / 1 , 


r-1 

^ n cos n^ = 
n=i 


r COB (r — 1)0 — (r — 1) cos r — 1 
2(1 — cos <fi) 
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If = 0, a numerical series is obtained 


= rjr --H 
2 


By use of these relationships in equation (7) it is found that, 
when k ± m is odd 




teik = 


2M. 


4nr sin — 
r 


1 


cos 


L 


f k + mW ]_ 
r 


— COB 


IlLrJSllI 


1 (8a) 


when k = m 


Pmk “ 


l 80 r 

8n Bin ^ 


and when k i m is even and k / m 


0mk " 0 


(8l>) 


(8c) 


For a symmetrical lift distribution 

= (211) 

V b \ b /r-ffl 


and 


°lk “Ir-k 


so that the summation for ou, needs to be made only from 
1 to r/2 
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r/2 

m=l 

where, when k ± m is odd 

Xmk = Pmk + Pr-m,k ® ^2) 

- l 80 

2nr sin -- 
r 

^ = 3mk ^ 

l8o 

itr(coB + 1) 

V r 


e ot _ cot 

ain Xk_t-ffibL sin 



(9) 


(10a) 


(lOh) 


when k = m 


Vk ~ ^mk 

l 80 r 

8jt sin ” 
r 

and when k ± m is even and k m 
\nk = ° 

For an anti symnetri cal lift distrihutlon 

/CjCN fcjc\ 

{-tl ' - 


(10c) 


(lOd) 
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and 


®lk “ “ 

In this case the summation 
to (I -l) since 

“Ik = 

vhere, when k ± m is odd 
^mk “ l^mk “ Pr-m,k 


for <Xijj needs to be made only from 1 
= 0; then 



l8o 

2itr 


sin2 -(^ 


Bin 


2 (k — m)?t 


(12a) 


when k = 


^mk “ Pmk 

l 80 r 


8it sin ^ 
r 


(12b) 


and when k i m is even and k m 


’'mk = 0 


(12c) 


Multipliers can thus be calculated so that the induced angle 
may be readily obtained by multiplying the known values of 

D 

by the appropriate multipliers and adding the resulting products. 
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The multipliers are independent of the aspect ratio eind taper ratio 
of the wing. Table e I and II present values of and 

and respectively, for r = 20. Similar tables for 

and are given in references 7 and 8, respectively, but 

no derivation is given therein. Tables for 

are given in reference h for valuee of r = 8^ l6, and 32, An 
inspection of totles I and II shows that positive valv.es occur on3y 
on the diagonal from upper left to lower right and that almost half 
of the values ai^e equal to ^ero. The multipliers and 

may he used with either nonlinear or linear section lift data whereas 
the multipliers for may be used only with linear section 

lift data. 

The method of determining the lift distribution becomes one 
of successive approximations. For a given geometric angle of attack, 
a dietrlbutlcn of c-^ is aaeumed from which the load distri— 

butlon is obtained. The induced angle of attack is then 

determined by equation (6), (9), or (11) through the use of the 
appropriate inultlpliers and subtracted from the geometric angle 
of attack to give the effectivo angle of attack at each spanvise 
station. From section data for the appropriate airfoil section and 
local FeyiiOlds number, values of ai’e road which correspond 

to the effective angle of attack of each section. If these values 
of c^ do not ag.'ee with those originally assumed, a second 
8SB\U£iptlon is made for and the process is repeated. Itirther 

assumptions cire made until the assumed values of are in agreement 

with those obtained from the section data. 


Wing Characteristics 

Once the lift distribution of a wing has been determined, the 
main part of the pi*oblem of calculating the wing characteristics 
is completed. The induced-drag and induced— yawing-moment coefficients 
are entirely dependent upon the lift distribution and it is assumed 
that the section profile-drag and pitching-moment coefficients are 
the same functions of the lift coefficient at each section of the 
wing as these determined in two-dimensional tests. 


The calculation of each of the wing coefficients involves a 
spanwise integration of the distribution of a particular . 

function I'hls Integration can be performed numerically 


through the use of additional sets of multipliers which are found 
in the following manner. 
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= ^ A, 
2 1 


/ 2v \ 

Since zhe values of A * yy are determined at the points 
can be found by hamonlc analysis as in equation (5) 



A 


1 


r-1 



m=l 


sin 


w 

r 


Therefore 



r-1 

\ 

= 2 \ 
m=l 


\ t /m 


•m 


Cl3a) 


where 
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If the distribution is sjTmnetriceLl,, f( and 

\ b \b /jv_3i 


where 


L 


2 y r(s^) 

0 h 




m 


4e18 


(13b) 


hms " 


im 


(^f\) 


^ms 


(-1) 


The moment of the dlst j'lbutlon f(~j can bo found in a similar 
manner. \b / 


4'^ 


£ <?) (?) <?) 


f}^ 




sin nd\ sin 6 cos 6 de 
/ 


= f 

4 2 


r-1 



(14a) 



NACA TN No . 1269 


15 


where 


= -5- 8ln 
m 8r r 


If the distribution is antisymmetrioal, f{ = — f(~) 

\b /jh \b /r— 12 


E- 3 . 


"-1 


m- 1 


(I4b) 


where 




Values of Cu,. '^ma are given in table III for r = 20. 

V^ing lift coef i ;1 dent The ving lift coefficient is obtained 
by means of a span-rflse integration of the lift distribution. 



If the lift distribution is esj-mmetrical 


r—1 



( 15 a) 
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If the lift dietrlbutlon is sj'iamctrlcal 


Cl = A 



^itia 


(15h) 


Ind uce d -drap coefficient The section induced-drag coefficient 

is equal to the product of the sect! on lift coefficient and the 
induced angle cf attack in radians. 



180 


The vring induced-drag coefficient is obtained by moans of a sparwiBe 
integra-vion of the section induced-drag coefficient multiplied by 
the local chord; 






180 ^ 


n-^ 


A 


2il nil 

b l 80 Vb*7 


For asymmetrical lift distributions 


-D; 


r-1 

rtA \ /CiC 

l83 2_ V b 

m=l 


h-ym 


^m 


(16a) 


For 8ymmetr.ical lift dictributions 
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P rof ile -dr 3 g_ coeffi ci e nt .- The esction profile-drag coefficient 

can le obtained from section data for the appropriate airfcil section 
and local ‘'teynolda number. For each spanvise station the profile- 
drag coefficient is read at tl;e aection lift coefficient proviously 
determined . The vring profile— drag coe'f'f icient is then obtained 
by means of a sparviae integration of the section profile -drag 
coefficient multiplied ly the local chord; 


(1V2 



1 

Q 

V.J 


I 



^do ^ 






For asyiTimetrical lift distributions 

r — 1 



or for syiiJiretrical lift distributions 



(17a) 


(iTt) 


P itchi ng-Eioment ^ofTlciert The eecticn pltching-moment 

coefficient a\>out its quarter-chord point can "be obtained from 
section data for the appropriate airfoil section and local 
Peynolds rju.niber'. For each spanwise station the pitching-moment 
coefficient is read at the section lift coefficient previously 
determined and then transferred to the wing reference point by the 
equation 
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“ ^o/k c 


COS 


(Og - C^ ) + Bin (ttg - tti)"] 

J 

C\ Bln (ag - ai) - cos (og - 0^)1 


( 18 ) 


vhere x and z are raeaaured froia the wing reference point to the 
quaj^ter— chord point of the section under coiiSideration and upvrard 
and haokward forces and dietencos aro taken as positive. The 
section pitclilng-Hncment coefficient about its aerodynamic center may 
be used instead of in which caeQ x and z aro meas’ired 


■■c 


to the section aerodynamic center. The terra sin (as - a^) may 

usually bo neglected . The wing pltching-Hraoment coef^’icient is obtained 
by the epanwlee integration 


Sc- 


,qV2 


-b/2 


Cm dy 


f 


1 Vcc*^ 

) \u ) 


For a8ymifletrical lift distrihutions 


r--! . 2>, 

^ V c c » L 

3 V C c /n 


( 19 a) 


For symmetrical lift distributions 




( 19 b) 


Eo lling-m om ent coeffici ent .- The rclling-moment coefficient 
is obtained by means of a spanwise integration 
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Cjcy dy 
Sb '' 

lJ-b/2 


_A ^ 

4 t b \b/ 

1-1 




For an anti symmetrical lift distribution 


2 . 


(¥). 


Induced-yawing-moment coefficient The induced-yawing- 

mcment coefficient is due to the moment of the induced— drag 
distribution 




Sb 1 l80 

U-b/2 


y 4y 


A £ic rto^i 2y 


..I..— . —■.I ■ 

b IBo b 


r— 1 

jtA / ^l9 


l8o ^ — \ ^ 

m=l 


H r. 
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The induced— yawing-Qoment coefficient for an antieymnietrical lift 
distribution is equal to zero and has little meaning inasmuch as 
the lift coefficient is also zero. The induced-yawing-nnoment 
coefficient is a f'unction of the lift and rolling-moment coefficients 
and must be found for asymmetrical lift distributions. 

Profile— yavlngHJOinent coefficient .— The profile— yawing- 
moment coefficient is due to the moment of the profile— drag 
distribution. 



n V2 
Sb 

U-b/2 


Cd^cy dy 



APPLICATION OF METHOD USING NONLINEAR SECTION LIFT DATA 
FOR SYMMETRICAL LIFT DISTRIBUTIONS 


The method desci-ibed is applied heroin to a wing, the geometric 
characteristics of which are given in table IV. Only symmetrical 
lift distributions are considered hereinafter inasmuch as these 
are believed to be sufficient for illustrating the method of 
cal dilation. The lift, profile-drag, and pltching-mcanent coefficients 
for the various wing sections along the span were derived from 
unpublished airfoil data obtained in the Langley two-dimensional 
low— turbulence pressure tunnel. The original airfoil data were 
cross-plotted against Reynolds number and thiolmoss ratio inasmuch 
as both varied along the span of the wing. Sample curves are given 
in figures 1 and 2 . From these plots the section characteristics 
at the various spanwise stations were determined and plotted in the 
conventional manner. (See fig. 3 *} The edge— velocity factor E, 
derived in reference 9 For elliptic wing, has been applied to 
the section angle of attack for each value of section lift coefficient 
as follows: 
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a„ = E 




- tt7 


+ 0.7 


Lift D1 str Ibvit ion 


Computation of the lift distribution at eJi angle of attack 
of 3° is shown in table V. This table is designed to be used where 
the multiplication is done by means of a slide rule or simple 
calculating machine. Wiere calculating machines capable of performing 
accumulative multiplication are available, the spaces for the 
individual products in columns (6) to ( 15 ) inay be omitted and the 
table made smaller. (See tables VII and VIII.) The mechanics of 
computing are explained in the table; however, the method for 
approximeting the lift coefficient distribution requires some 
explanation. The Initially assur.ied lift— coefficient distribution 
(column ( 3 ) of first division) can be taken as the distribution 
given by the geometric angles of attack but it is be.st determined 
by some simple method which will give a close approximation to the 
actual distribution. The initial distribution given in table V 
was approximated by 



where coefficient read from the section curves 

for the geometric angles of attack. This equation weights the lift 
distribution according to the average of the chord distribution of 
the wing under consideration and that of an elliptical wing of the 
same aspect ratio and span, tdien the lift distributions at several 
angles of attack are to be computed and after they have been obtained 
for two angles, the initial assmned cj distribution for subsequent 
angles can be more accurately estimated In the following manner; 
Values of downwash angle are first estimated by extrapolating from 
values for the preceding wing angles, and then, for the resulting 
effective angles of attack, the lift coefficients are read from the 
section curves. 

The lift coefficients in column (l3) of table V, read from 
section lift curves for the effective angles of attack, will usually 
not check the assumed values for the first approximation. In 
order to select assumed values for subsequent approximations, 
the following simple method has been found to yield satisfactory 
results. An Incremental value of lift coefficient Is obtained 

according to the relation (numbers in parenthesis are columns 
in table V) : 
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_(l8) — ( 3 )jnt- 3 _ + 3 . (^^) ~ ( 3 )_ m Tc^^) ~ ( 3 ) hh-1 

K 


where K has the following values at the spanwise stations 



K 

b 


0 to 0.8910 

8 to 10 

.9511 

11 to 13 

.9877 

l 4 to 16 


and (l8) - ( 3 )Jm is the difference between the check and assumed 
values' for the mth spanwise station. The incremental values so 
determined are added to the assiuned values in order to obtain new 
assumed values to be used in the next approximation. This method 
has been found in practice to make the check and assumed values 
converge in about three approximations if the first approximation 
is not too much in error. 


Wing Coefficients 

Computations of the wing lift, profile-drag, induced-drag, 
and pitching-moment coefficients are shown in table VI. Since the 
lateral axis through the wing reference point contains the quarter- 
chord points of each section, the x and z distances in 
equation (l8) are zero, and the pitching-moment coefficient of 
the wing is determined solely by the values of Cm /■ . 


APPLICATION OF METHOD USING LINEAR SECTION LIFT DATA 
FOR SYMMETRICAL LIFT DISTRIBUTIONS 


Although the method described herein was developed particularly 
for use with nonlinear section lift data, it is readily adaptable 
for use with linear section lift data with a resulting reduction in 
computing time as compared with most existing methods. When the 
section lift cui-ves can be assumed linear, it is usually convenient 
to divide any symmetrical lift distribution (as in reference 10) 
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into tvo parts — the additional lift distribution due to angle of 
attack chaagos and the basic lift distribution due to aerodynairic 
twist. The calculation of these lift distributions is illtinbrated 
in tables VII to X for the wingjtJtie gecsaetrJc characteristics of 
which were given In table IV. 

It shovild be noted that tables VII and VIII are essentially 
the same as table V but are designed primarily for use with 
calculating raachines capable of performing accumulative multi- 
plication. If such machines are not available, these tables r.iay 
be coii8bi"acted similar to table V to allow spaces for writing the 
Individual products. 


Lift Characteristics 


Tvro lift distributions are required for the determination of 
the additional and bae^c lift distributions. The first one is 
obtained in table VII for a. conetamt angle of attack Og (c ' C) 

and the second one in table VIII for the angle of attack distribution 

c c 

due to the aerodynamic twist (aag = O) . The check values of — ^ 
(coluitin (18)) are obtained by multlplj’-ing the effective angle 

of attack Oq . by final apprcocimations are entered in 

table rx as and 

V b A- s \ ^ Jf 

The — \ distribution is the additional, lift distribution 

'' 

corresponding to a wing lift coefficient Cp,/ . determined in 

vccas ) 

table IX through the use of the multipliers It is usually 

convenient to use the additional lift distribution - corresponding 

to a wing lift coefficient of unity. This distribution is found 
by dividing the values of (——) Ct , 

k ^ /(“as) 


/ct c \ 

The V -r- I dif^trihution is a comhination of the hasic lift 

\ b /(e^>) 

distribution and an additional lift distribution corresponding to 


a wing lift ccefficieut Cp 




(n*) 


also detormined in table IX. The 


basic lift distribution -i;~ is then determined by subtracting the 

^ c, _c /c,c' 


additional lift distribution — 


^al' 


Cl 


(H') 


from 


\h')‘ 
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Inaeinuch ao the wing lift curve is asBumed to he linear, it is 

defined By ita slOiPe and angle of attack for zero lift vhich are 

alao found in table IX. The maximim wing lift coefficient is 

estimated according to the method of reference 10 which is illustrated 

in figure 4. The maxinum lift coefficient is considered to be the 

wing lift coefficient at which some section of the wing becomes 

the first to I'each its maximum lift, that is, Ci, + Ct ci = 

^h ‘mI ^laex 

This value of Ct is moot conveniently determined by finding tie 

minimum, value of — ^2. along the span as illustrated in table IX. 

^iel 


Induced— Drag Coefficient 

The section Indi' end-drag coefficient is equal to the product 
of the section lift, cenff icient tmd the induced angle of attack 
in radians. The lift distribution for any wing lift coefficient is 


S'jJi = 
b b 


Cl + 



(23) 


The coiTeapondlnf' induced rngle of attack distribution may be 
vribten as 


Cl + C^,-^ (24) 

The values of and are determined in table X in the 

same manner as — pnd — ii- in table IX. The induced-drag 
b b 

distribution is therefore 

*^dj° _ aq 

~ ~ 'b~ 57-3 


^di ^ 
b 


"dial" 2 

^ Cl^ + 


alb 


Cl + 


cdi^c 


(25) 
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where 


t b 57.3 


(26) 




a. 


J.i^ 

57.3 



a.f ^ 

57.3 


(27) 


and 





(26) 


The caTcu?.ab:' o:"! of each of ifieoe iiiduced-drag di3tri Out ions Is 
illustr*abv d 'd. tahlc X together with the nurnorioal integration of 
each d-strihivtion to obtain the wing indue ed-drcog coefficient, 

?roflle-frag and Pitching— MoMent Coefficients 

The rrof ile-drag and pitching-iTonent coefficients for the 
wing dopond rirect]p; TUjon the section data and therefore their 
ca.lcula;l' icn is the sai.io whether linear or nonlinear section lift 
data are used. For tae linear case the section lift coefficient is 


— c 1 


Cl 


for any wing co-^fflcient Cj^. By use of this 

j»rof;llG-irttg and pitching— sioment coeff Icienta 
table Yl. 


value for 
are found 


c 2 the 
as in 


DISCUGSTON 


The char<j.ctorl3tics of three wings w?.th syinrtetrical lift 
distributions have been calculated by use of both nonlinear and linear 
section lift data and are presented In figure 5 together with 
experimental results. The bo data were taken from i-eferencs 11 . 

Ti^e lift curves calculated by use of nonlinear section lift data 
ar-e in close agreeTPient with the experlmeiital results over the 
entire range of lift coefficients- whereas those calculated by use 
of linear section lift date are in agreement only over the linear 
portions of the ciurves as would he expected. 
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It nr.iBt be remembered that the methods presented are subject 
to the limitations of lifting-line theory upon which the methods 
are based; therefore, the close agreement shown in figure 5 should 
not be cTccected for wings of Icni aspect i-atio or large sweep. 

The use of the edge— velocity factor more or lees compensates for 
some of the effects o:^ aopect ratio and, in fact, appears to over 
con 5 )ensate at the larger values of aspect ratio as shewn in figure 5 . 

Additional comparisons of calculated and experimental data 
are given in reference 11 for wings with syLimetrical lift distributions, 
but very little compjarable data are available for wings vrith 
asymmetrical lift distributions. Such dt-.ta are very desirable 
in order to determine the roliability with which calculated data 
may be used to predict experimental wing choracterietics. 


Langley Memorial Aeronautical Laboi’etory 

Kationai- Advisory Commlttoe for Aeronautics 
Langley Field, Va. Decembei' 20, 19hS 
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Values of k at top to be used with values of m at left side. NATIONAL ADVISORY 

COMMITTEE FOfl AERONAUTICS 

Values of k at bottoin to be used with values of m at right side. 
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TABLE III.- WING-COEPFICIENT MULTIPLIERS 


¥ 

m 


^3 



-0.9877 

19 

0.01229 


-0.00607 


-.9511 

18 

.021+27 


-.01151+ 


-.8910 

17 

.05066 ^ 


-.01589 


-.8090 

16 

.01+616 


-.01867 


-.7071 

15 

.05551+ 


-.01961+ 


-.5878 

14 

.06551+ 


-.01867 


"•1+5U0 

15 

.06998 


-.01589 


-.3090 

12 

.071+70 


-.01151+ 


-.1561+ 

11 

.07757 


-.00607 


0 I 

1 

10 

.07851+ 

0.07851+ 

0 

0 

.156[|. 

9 

.07757 

.15515 

.00607 

O.OI21I+ 



.5090 

8 

.071+70 

. 11+959 

.01151+ 

.02508 

.1+81+0 ! 

1 j 

7 

.06998 

.15996 

.01589 

.05177 

.5878 

6 

.06551+ 

.12708 

.01867 

.05755 

.7071 

5 

1 

.05551+ 

.11107 

.01961+ 

.05927 

.8090 

1 

1+ ! 

1 

.OI+616 

.09253 

.01867 

.05735 

.8910 

5 

-^.05666 

.07151 1 

.01589 

.05177 

.9511 

2 

/ .021+27 ' 

.01+851+ 

.01151+ 

.02508 

.9877 

1 i 
i 

.01229 

.021+57 

.00607 

.OI21I+ 
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TABLE I 
First opproximotion 


- CALCULATION 


DISTRIBUTION 


II) 

(2) 

13) 

(4) 

(5) 

( 6) 

(7) 

(8) 

( 9} 

( 10) 

(M) 

( 12) 

(13) 

( 14) 

( 15) 

(16) 

(IT) 

(18) 


a 

<a 

c 

b 






^mh * 

column (5) 




O'i 

«e 



(oe,* el 

louurncd) 

(Tobte U) 

0x14) 

0 

.1564 

.3090 

.4540 

.5878 

.7071 

.8090 

8910 

.9511 

.98 77 

S(6) 

k>(i9» 

I2)-(I6) 

(cl'eciO 






143.239 

-50.533 

0 

-6.950 

0 

-2 865 

0 

-1.804 

0 

-1.466 





^.0( 

Q.Sl’ 

0.lii2< 

0.075 

i 10. SO 

-li.29 

0 

-.*>1 

0 

-.21 

0 

-.15 

0 

-.11 

1.88 

1.12 

0 .Ubii 






-115.624 

145.025 

-67298 

0 

-10.158 

0 

-4 840^ 

p 

-3 394 

0 




0. 1 5 6 4 


»S17 

.129« 

.067C 

-7.7S 

9.72 

-U.si 

0 

-.63 

0 

-.52 

0 

-.25 

0 

.98 

1.78 

.S5] 






0 

-64.802 

150 6 II 

-67.157 

0 

-9 916 

0 

-4968 

0 

-3.768 




.3090 

ZM 



. 060 s 

Q 

-5.9S 

9.17 

-I 1.09 

0 

-.60 

0 

-.50 


-.25 

.90 

1.S7 







-12.384 

0 

-62.917 

160.76 1 

-72.4 72 

0 

-10 926 

0 

-5.812 

0 




.4540 

2.;i 

.Sl<5 

.lOiiC 

.05liC 

-.67 

0 

- 5 .I 1 O 

8.60 

-5.91 

0 

-.S9 

0 

-.51 

_Q 

^77 

1.56 

.•ii 


! 




0 

-8 320 

0 

-65.803 

177.054 

-82083 

0 

-13.134 

0 

-7.713 




.5878 

>1-73 


i -Q?2! 

.0U6! 

0 J 

_--12 , 

0 

- 3.05 

6.20 

- 5.80 

0 

-.61 

0 

-.56 



. 5 a 






-4.05! ^ 

0 

-7.372 

0 

-71.743 

202.571 

-97965 

0 

-17 388 

0 




.7071 


.U77 

. 062 > 

.059! 

- .16 ^ 

0 

-.25 . 

0 j 

- 2.82 

7.96 

^S5. . 

0 

-.68 

0 

-.65 

i *63 

, . 4 -^ 






0 

-2.880 ^ 

0 

-7.208 

0 

-81.434 

243 694 

^25.537 

0 

-26.635 


1 

[ 

.8090 

M 

•U3C 

.075' 

.051^ 

0 

>.09 

0 

-.23 

0 

- 2.57 

7.70 

-5-97 

0 

..aiL 

.2s 

_.-4^ 






-1638 


-2 371 _ 

0 

-7.370 

0 

-96 962 

315.512 

-I80.52ff 

q 




.8910 

_ J2 

j6c 

,0663 

__ ^0239 

-.(4 

0 

-.06 

0 

-.18 

0 

- 2.32 


- L .51 

0 

*42 

-.10 

•4^i 






0 

-1 06 2 

0 

-2 016 

6 

-7 599 

0 

-122.886 

463J533 

-329976 




.9511 

-.10 

. 28 J 

. 061 ! 

.0172 

0 

-.02 , 

0 

. -.03 

0 

- -±li . - 

0 

- 2.11 

7.97 1 

-5.66 

•77 

-.87 







-0 459 

0 

-0,620 

0 

-1.491 

0 

-7039 

0 

-167045 

915 651 




.9877 

I - .39 


1 «Qji51 


0 

p ^ 

-.01 

0 

-.01 

0 

-07 

0 

- 1-67_J 

.9-l6_ J 

Ri 

- 2.33 _ 

. 16 F 


Second approximation 


2, hi .5^4 .116L .0610 
.«^I7 -lOUQ .0S58 


-58.533 0 

-6950 

-4.17 , 9 .. ! 

-.li9 

145 025 -67.298 

0 

9.69 -li.SO 

0 

-64 802 I506H 

-67167 

-3.95 9.19"^ 

♦4.10 

0 -62.917 

1 6076j 

0 -5.38 

8.6s 

-8 320 0 

-65803 


-10 158 *^0 



-3394' 

tap , 

0 

1.61 

1.59 

-.23 

0 

0 

-3 76 8 

1.07 


0 i 

-.2~5 

•?5 

1-52 

-5. 812 1 

0 




..232LQ613-. 4^511 








611 -6 7.15 7 O 

9.16 -L.oe 0 -.60 ^1 

62.917 160 76lp727472 0 j 

5.38 8.63 -3.89 p 1 

O -65 803 177.054-82.683 


-.10 , 

l.SS 

l.liS 

•ii 97 

0 




0 I 

1,12 

1.6L 

.Slfi 

-3.768 




-.23 

.91 

1.S6 

.S 21 

q 

0 


1.39 

.“jl? 

-7 713 

.60 , 

1.15 

.SOO 

' 0^ _ 




0 

.59 

..69 

.479 

- 26535 ^ 

-. 62 ^ 

-ilS , 

.442 


H67045 915 


Numbers oppearing in parentheses denote column number. 
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NACA TN No. 1269 


Fig. 1 







NACA TN No. 1269 




Fig. 3b 
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NACA TN No. 1269 


Fig. 3c 



(c) Pltchlng-moment . 


Fig. 4 


NACA TN No. 1269 



Figure Estimation of for example wing. estimated to be 1.37-) 


Fig. 5a 



(a) A= 8.0U, R= 4»320,000, root section NACA tip section NACA 

Pigure 3«» Experimental and calculated characteristics of three wings of taper ratio i 
and NACA l^i^-series airfoil sections. 


Fig. 5b 


NACA TN No. 1269 



(b) A= 10 . 05 , R= 3 » 490*000, root section WACA 4420, tip section KACA 4412. 
Figure 5»- Continued. 




Fig. 5c 



12.06, R= 2,870,000, root section NACA 14^4214., tip eectlon MACA U|12. 





